Transport in a Levy ratchet: Group velocity and distribution spread 
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We consider the motion of an overdamped particle in a periodic potential lacking spatial symmetry 
under the influence of symmetric Levy noise, being a minimal setup for a "Levy ratchet." Due 
to the non-thermal character of the Levy noise, the particle exhibits a motion with a preferred 
direction even in the absence of whatever additional time-dependent forces. The examination of 
the Levy ratchet has to be based on the characteristics of directionality which are different from 
typically used measures like mean current and the dispersion of particles' positions, since these get 
inappropriate when the moments of the noise diverge. To overcome this problem, we discuss robust 
measures of directionality of transport like the position of the median of the particles displacements' 
distribution characterizing the group velocity, and the interquantile distance giving the measure of 
the distributions' width. Moreover, we analyze the behavior of splitting probabilities for leaving 
an interval of a given length unveiling qualitative differences between the noises with Levy indices 
below and above unity. Finally, we inspect the problem of the first escape from an interval of given 
length revealing independence of exit times on the structure of the potential. 

PACS numbers: 05.40.Fb, 05.10.Gg, 02.50.-r, 02.50.Ey, 



I. INTRODUCTION 

Motion of particles in an external periodic potential 
with broken spatial symmetry and under influence of ex- 
ternal forces with zero mean can result in occurrence of 
the persistent, directed current [l], 0, S 0, III • The occur- 
rence of the directed transport is a consequence of the 
violation of the detailed balance in systems acting away 
from equilibrium. Typically, one assumes that the overall 
force acting on the particle is a superposition of the Gaus- 
sian thermal noise with another periodic or stochastic 
force [i} . In those cases all moments of the distribution 
of the noisy force do exist, and the existence of moments 
of the distribution of the particle's velocity is guaranteed. 
There are at least two reasons to consider the directed 
motion of a particle in a periodic ratchet potential under 
the influence of a heavy-tailed noise causing anomalously 
large particle displacements. One of them is pursuing the 
line of investigations of rectification of non-thermal noises 
and generalizing the corresponding considerations to the 
heavy-tailed cases which were found to be quite abun- 
dant under non-equilibrium conditions (y, [7|, |8|, l2( . An- 
other one, connected with the former, is seeking for the 
way to characterize the ensuing directed motion in the 
case when the corresponding statistical moments are ab- 
sent, so that neither the dispersion nor even the mean of 
the corresponding displacements exist. In consequence, 
in such cases the standard characteristics of motion like 
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mean velocity or Peclet number become inapplicable. 

The interplay of deterministic dynamics and pertur- 
bative Levy-type noises have been addressed in litera- 
ture in various scenarios including several noise-induced 
effects like resonant activation [bCA Jill , stochastic res- 
onance [12], dynamical hysteresis [l2|, LL3(, studies of 
decay/relaxation properties of the probability densities 
[l4 Il5| , escape from bounded interva ls [la . I17j , the clas- 
sical barrier crossing proble m 118. [l9l . 120 ] or examination 
of stationary states [la, [2J, [22|, [23| . However, very few 
examples [24| tackle the problem of Levy noise driven 
dynamics of periodic systems. 

In the present work we study the behavior of a particle 
in the periodic potential with a broken spatial symme- 
try subjected to the action of a symmetric Levy stable 
noise. The a-stable noise can be expected to occur in 
systems being out of thermal equilibrium, where detailed 
balance conditions is violated. In such realms one can an- 
ticipate noise induced current in a static potential with 
broken spatial symmetry. Consequently, what we ad- 
dress here as the Levy ratchet resembles so called ther- 
mal ratchet 0] with time-dependent temperature varia- 
tions. Heavy tails of distribution of additive stochastic 
increments lead however to considerable peculiarities of 
such motion, due to power-law character of probability 
densities of the particle's displacements, which lack the 
dispersion and may also lack the mean. Therefore, the 
examination of the current defined as the time deriva- 
tive of the mean position may not be adequate for the 
Levy ratchet, and one has to look for other quantities 
which characterize the motion caused by the interplay 
of the noise and a potential. Instead of examination of 
the current and of the effective diffusion coefficient we 



discuss several robust probabilistic characteristics of the 
particle's displacement based on its cumulative distribu- 
tion: the behavior of median, which allows to introduce 
the group velocity, the growth of interquantile distances 
and the behavior of splitting probabilities (as defined in 

sec. hud). 

As it will be discussed futher on, symmetric stable 
noise acting in a presence of a potential with broken 
spatial symmetry can induce directed motion of a test 
particle in an overdamped system. As long as the mean 
value of the perturbative noise exists and equals zero the 
displacement of the particle also possesses a mean. The 
overall motion can be then characterized by the temporal 
changes in the mean position, or, in the other words, by 
a mean velocity. In the case of an ensemble of particles 
mean velocity corresponds to the overall current, which 
is an adequate measure of the directionality of motion 
[J]. This situation is realized for a-stable noises with 
1 < a< 2. 

If the noise, however, is distributed according to a 
Levy-stable law with < a < 1, neither the mean of 
the noise nor the mean of the overall displacement exist, 
so that one has to look for other quantities characteriz- 
ing the directed motion. The discussion of these quanti- 
ties, parallelly to investigation of the properties of the di- 
rected transport, is the main topic of the present article. 
Its structure is as follows. The next section (Section |TT) 
presents the model under consideration. The obtained 
numerical results are included in Section UTIl which is de- 
voted to various characteristics of the considered ratch- 
eting device. The cumulative distribution of the over- 
all displacement probability is examined in Section [ill Al 
whereas its main features like median allowing for a def- 
inition of a group velocity of the probability packet are 
studied in Section IIIIBl We discuss also interquantile 
distance (Section IIIIC|) . as well as splitting probability 
(Section IIII D[) that defines the fraction of current going 
to the left versus right direction when leaving an interval 
of finite length. Finally, (Sec. lIIIE) we discuss properties 
of the first escape time distribution from a box of given 
width. The paper is closed with summary and concluding 
remarks (Section IIV[) . 



II. MODEL 



Initially, at t — 0, a test particle is located at x = 0. 




FIG. 1: The potential, V(x), with broken spatial symmetry 
given by Eq. ((2) used for examination of the studied Levy 
ratchet. 



In the most general case Levy distributions correspond 
to a 4-parameters family of the probability density func- 
tions [25|, |2(| |27| . Here we restrict studies to symmetric 
stable noises, L a ((; a, //), which are characterized by the 
Fourier-transform 4>(k) = f^° e lfcl »L a (C; ^M^C of the 
probability density given by [25|, |26J, |23| 



4>{k) = exp[ikfi-a a \k\ a ] 



(3) 



The parameter a (where a € (0, 2]) is the stability index 
of the distribution describing (for a < 2) its asymptotic 
"fat" tail characteristics yielding L a (£; a, fi) ~ |£|~( 1+Q ) 
for large (. The parameter fi denotes the location pa- 
rameter representing position of the modal value. The 
Gaussian distribution, see Eq. ([3), corresponds to a spe- 
cial case of a Levy law for a = 2, with /i interpreted now 
as a mean and a as the dispersion of the distribution. 
Additionally, we assume \i = and a — {0.25, 0.5, 1}, i.e. 
we consider symmetric, strictly stable noises only. It is 
worthy mentioning that nonzero asymmetry parameter, 
which is not considered here, introduces a preferred di- 
rection of noise pulses and consequently may lead to a 
preferred direction of motion [23l . |24| |. 



In what follows we consider the system described by 
the overdamped Langevin equation 



I— V W+ ««), 



(i) 



where £(£) stands for a Levy stable noise, i.e. the noise 
for which increments are distributed according to a sta- 
ble probability density. Finally, V(x) is an external static 
ratchet potential, see Fig. [T] with broken spatial symme- 
try 
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III. RESULTS 

For symmetric Levy noises, Eq. ((T) is equivalent to the 
following fractional differential Fokker-Planck equation 

nnn 

0P fo*> = °-V'(x)P(x t) + o* * P W (4) 

dt a/ lJl,, + ff d\x\ a ' ( ' 

where a and a characterize the noise. The frac- 
tional (Riesz-Weyl) derivative is interpreted in the sense 
of the Fourier transform HH HI HH J^f{x) = 



aRF 



■II 



dk, 



-ikx 



\k\ a f(k). Nevertheless, due to possible in- 



stabilities [32] of numerical approximations to Eq. (j4|) [33] 
in the following studies we use approach based solely on 
the Langevin equation {J). 

The results presented below were obtained by the nu- 
merical integration of Eq. (JXJ) with Levy stable noises 
using the methods discussed in Refs. jig, [2QJ, 12s l27j ]. 
The numerical integration of Eq. JT]) was performed with 
the time step of At = 10~ 3 . The number of realizations 
varied from N — 10 5 to N = 10 6 leading in all cases to 
consistent results. 

Initially, a test particle is located at x = 0. In course of 
time the distribution of the particle's positions x evolves 
due to the presence of the stochastic force and of the 
deterministic potential in Eq. (fTJ) . As we proceed to 
show, the ensuing coarse-grained probability distribution 
of a; is a symmetric distribution of the power-law type 
(Sec. IIII A[) . The overall motion of the probability den- 
sity is characterized by changes in the location of its me- 
dian (Sec. IIIIBp . The width of the probability density 
grows with time. The growing width of the distribu- 
tion is characterized by growing interquantile distance 
(Sec. IIII C|) . Still another characteristic of the motion 
can be addressed by analyzing behavior of splitting prob- 
abilities: A particle initially located at x = can escape 
from an interval of given width L centered at the ini- 
tial position to the left or to the right. The statistics of 
the corresponding escape events (splitting probabilities, 
Sec. IIII D|) provides another measure of the directional- 
ity of motion. Finally, mean first escape time provides a 
characteristic of the time scales involved in such escapes 

(Sec. mm). 



A. Cumulative distribution 

We first discuss the cumulative distribution func- 
tion (CDF) of the particle's displacement under the 
non-equilibrium additive noise. The CDF is given by 
F(x,t) = f* P(x',t)dx' with P(x,t) being the corre- 
sponding probability density function (PDF), which can 
be interpreted in terms of particle concentration. Simula- 
tions show, that in course of time the corresponding PDF 
attains a symmetric form, see Fig. [5] The density P(x, t) 
is then found to move in one direction and to broaden 
in time. These effects can be seen by inspection of the 
overall position of the particles characterized by the me- 
dian of the distribution q , 5 (cf. Fig. [5]) and its width 
defined by the interquantile distance, e.g. qo.s{t) — 90.2 (£) 
(cf. Fig. 0). The quantiles q m (0 < m < 1) are defined 
according to the relation F(q m ,t) = m. 

For symmetric probability densities, i.e. densities such 
that P{qo.5 —x,t) = P(x + qo.s,t) the cumulative density 
F(x,t) fulfills 



F(x,t) = l-F(2q . 5 -x,t). 



(5) 



spect to median (90.5) if the space scales much larger 
than the period of the potential are considered, see Fig.Q] 
and [2] 
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For the Levy ratchet system, the constructed CDFs prove 
that PDF recorded at given time are symmetric with re- 



FIG. 2: (Color online) Cumulative distributions F(x,t — 100) 
and 1 — F(2qo.5 — x,t — 100) for a = 0.7 (top left panel), 
a = 0.9 (top right panel), a — 1.5 (bottom left panel), a — 1.7 
(bottom right panel) with a — 1. The perfect agreement 
between both curves indicates the symmetric shape of corre- 
sponding PDFs with respect to medians. The same agreement 
is observed for smaller values of the noise intensity a. Note 
the large difference in scales for different values of a. 



The distributions of the particles' displacements in 
Fig. [2] exhibit power-law tails. To detect them it is 
enough to consider the asymptotic behavior of the 'sur- 
vival probability' G(x,t) = 1 — F(x,t) for large x, see 
Fig. [3] clearly displaying this power-law behavior. The 
simplest hypothesis here would be that the correspond- 
ing exponent is the same as the one of the noise; if £(£) 
has a stable density characterized by exponent a, i.e. 
P((,t) ex e~( 1+Q ) for large £> the 'survival probability' 
should behave like G(x,t) ex x 1 with 7 = a. This kind 
of behavior is indeed observed for the values of stability 
index a smaller or around one, see Fig. 0] For larger 
a, G(x, t) still shows a power-law asymptotic, see Fig. [3[ 
however, the numerically obtained values of exponents 7 
differ from a, see Fig. [4] This difference is pertinent to 
a very far tail of the distribution, since the behavior of 
the interquantile distances, say go. 2 {t) — q_o.a{t), unveils 
the behavior compatible with the 7 = a assumption for 
all a, see Sec. IIII CI 

Inspection of 'survival probability' G(x,t), see Fig. 02 
suggests that for small value of a the particle practically 
does not feel the potential. Therefore, the particle be- 
haves like a free particle and the exponent characterizing 
'survival probability' is the same as the one characteriz- 
ing statistical properties of the noise. This observation 
is in accordance with predictions of the continuous time 
random walks theory in situations when the mean wait- 
ing time for a next jump is finite [34 , l35| . For increas- 
ing values of a, the motion gets to be more sensitive to 
the structure of the potential. Consequently, the form of 



G(x, t) in the intermediate range of x diverges from the 
pure Levy distribution. The difference between 7 and a 
is probably explained by extremely slow convergence of 
the tail of the corresponding distribution to its asymp- 
totic form. 

In order to fit the tail asymptotics properly, it is nec- 
essary to have large statistics corresponding to extreme 
events which should cover at least a couple of orders of 
magnitude. Such a coverage is rather easily reached for 
values of a ^ 1.1. However, for large values of the sta- 
bility index a (a > 1.1), the statistics of 'survival prob- 
ability' G(x, i) for very large values of x is quite poor: 
the exploration of large x dependence requires simulation 
times which are beyond our possibilities. The values of 7 
obtained in simulations did not indicate changes for cho- 
sen simulation times, however the fitted slope was found 
to be sensitive to the choice of the threshold from which 
a power-law is fitted. 




FIG. 3: (Color online) 'Survival probability' G(x, t = 100) = 
1 — F(x,t — 100) with a = 1. For large x, numerical simula- 
tions demonstrate power-law character of 'survival probabil- 
ities' G(x,t = 100). Nevertheless, the value of the exponent 
characterizing large x behavior is different from a when un- 
derlying noise is characterized by a^ 1.3, see Fig. [4] 



B. Median and group velocity 

The position of the median describes the overall mo- 
tion of the probability density to find a particle in the 
vicinity of x. In our case, the corresponding PDF coarse- 
grained over the period of the potential, is a symmetric, 
monomodal function, so that the temporal change of the 
median (which coincides with the maximum of the PDF) 
gives us the possibility to define the group velocity of the 
particles' packet as v = dqo.s/dt. Our numerical simu- 
lations show that the position of median go.s(i) changes 
linearly with time 
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FIG. 4: Value of the exponent 7 in 'survival probability' de- 
scribing the asymptotic behavior of its tails G(x,t) oc x 1 , see 
Fig- E] For q ^ 1.3, the exponent characterizing tails' asymp- 
totic differs from the index a of the underlying noise because 
the motion of the test particle starts to be more sensitive to 
the structure of the potential. Relative errors of the fit are 
smaller than 1% of value of exponents. Consequently, error 
bars are smaller than symbol size. 



90.5 (t) 



t + b, 



(6) 



which corresponds to the constant group velocity. Note 
that this group velocity can be defined even in the case 
when the mean velocity or current does not exist due to 
large fluctuations, and defines the behavior of the typi- 
cal displacement of particles with time. Fig. [5] presents 
location of the median as a function of time for different 
values of stability exponents a. The 'periodic-like' mod- 
ulation of the median position is caused by the periodic 
shape of the potential, see ordinate of Fig. Q] For the 
constant amplitude of the noise the group velocity de- 
creases with increase of the exponent a and disappears 
in the Gaussian case a = 2, see Fig. [6] in which case the 
noise satisfies the detailed balance condition. 

The 'survival probability' shows that asymptotic be- 
havior of the probability density of a particle position is 
not sensitive to the detailed structure of the potential. 
The presence of potential is visible for larger values of 
a when the exponent characterizing the tails of distribu- 
tions starts to differ from a, see Figs. [3] and [4j Different 
situation takes place for median which reflects the struc- 
ture of the potential. It is manifested in the periodic 
modulation of median position visible in Fig. [5j This pe- 
riodic structure has the same period like the potential. 
This suggests that tails of distributions can freely over- 
pass the potential maxima, while the motion of median is 
strongly affected by presence of maxima of the potential. 
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FIG. 5: (Color online) Time dependence of the location of 
median, qo.s(t) which describes the overall motion of the en- 
semble of particles injected into the Levy ratchet at t — 
(<r = 1). The linear growth of median with time allows to 
define the group velocity of the packet of particles. 
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FIG. 6: Value of the group velocity (coefficient v in Eq. (|6} 
characterizing slope of curves in Fig. [5} as a function of the 
stability index a. The lines are drawn to guide the eye only. 
The fastest overall motion of the probability density is ob- 
served for small values of stability index a for which large 
noise increments are more probable. The value of v van- 
ishes towards a = 2 corresponding to the equilibrium Gaus- 
sian noise for which persistent current disappears. Error bars 
represent errors of the fit. Error bars, except for those for 
a = 0.5, are smaller than symbol size. 



C. Interquantile distance 

Interquantile distance is a robust measure characteriz- 
ing the distribution's width even in cases when moments 
of the distribution do not exist. For the studied Levy 
ratchet this interquantile distance scales like t p indepen- 
dently of the value of the noise intensity a, i.e. 



where p = 1/a. Fig. [7] presents time dependence of 
interquantile distance qo.2 (t) — qo.s(t) along with fitted 
curves. Fig. [5] presents values of fitted exponents, to- 
gether with the predicted 1/a-behavior. The interquan- 
tile distance q m (t) — qi- m (t) scales in the same manner 
also for other values of to (0 < to < 1). Note that the 
periodic modulation which is typical for evolution of me- 
dian, see Fig. [5J can be not visible for interquantile dis- 
tance. Essentially, this modulation does not appear in 
the behavior of whatever quantiles q m (t) with to =/= 0.5 
and sufficiently large noise intensity a. The decrease of 
the noise intensity a, introduces periodic modulation also 
to interquantile distance which can be the most easily ob- 
served for large values of the stability index a (figure not 
shown). 




Qm(t) ~ qi-m{t) OC t p , 



(7) 



100 



FIG. 7: (Color online) Time dependence of the interquantile 
distance qo.2(t) — qo.s(t) with a — 1 together with fitted the- 
oretical curves. The decrease of scale parameter a makes the 
particle more sensitive to the shape of the potential and con- 
sequently introduces periodic- like modulation of interquantile 
distances (results not shown). 



D. Splitting probability 

Considering splitting probabilities gives us another 
characteristics of the directionality of motion. Splitting 
probability is a measure of the preferred direction to 
which the first escape from the system takes place. In this 
case we consider an interval of length L centered around 
the point where particles are introduced, and consider the 
probability for the particle to escape the interval through 
its left (right) boundary. Fig. [5] presents probability tt^ 
of the first escape to left as a function of the box half- 
width. The presented splitting probability indicates that 
the preferred direction of the first escape from the system 
is to the right, since for longer intervals ttl < 1/2. Fur- 
thermore, decrease of the scale parameter a makes this 
tendency stronger. 

Fig. [5] shows a clear difference between the cases 1 < 
a < 2 and a < 1. In the first case for the intervals 
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boundary. 

In the case a < 1 the width of the distribution 
grows faster than its maximum moves. Due to the well- 
prescribed initial position, the contribution of the overall 
motion is still perceptible at shorter times and thus for 
smaller boxes, and the difference between the probabili- 
ties to leave the box through its left and right boundary 
grows. At longer times the broadening of the distribu- 
tion prevails, typical displacement gets hardly relevant 
on the background of fluctuations, and the difference in 
splitting probabilities decays. However, the distribution 
still has a pronounced maximum, and its motion defines 
the group velocity. 



FIG. 8: Value of the exponent p in Eq. ((7} as a function of the 
stability index a. Line presents theoretical dependence of the 
exponent characterizing growth of the interquantile distance, 
i.e. 1/a. Error bars are smaller than symbol size. 
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FIG. 9: (Color online) Splitting probability, ttl, i.e. proba- 
bility of the first escape to the left as a function of the box 
half-width. Initially particle is located in the middle of the 
box of half-width L/2, see Eq. @ and Fig. [T] Various pan- 
els correspond to different values of the scale parameter a: 
a = 0.25 (left panel), a — 0.5 (middle panel) and a — 1 (right 
panel) . 



("boxes") much larger than the period of the potential 
the probability ttl decays monotonically with the box 
size, which means than more and more particles leave 
the box through its right boundary. In the second case, 
as exemplified by a — 0.9 the behavior of ttl is nonmo- 
motonical: it decays with L up to L/2 w 20 and then 
starts to grow again. 

This striking difference between the results for a > 1 
and ones for a < 1 gets clear when returning to the 
behavior of the interquantile distances. In the first case 
the shift of the distribution's median (coinciding with 
its maximum), going linearly in time, is asymptotically 
faster than the growth of the distribution's width, which 
allows us to discuss the position of the typical particle as 
a mean position and to introduce the mean current. In 
course of the time (necessary to leave the larger box) the 
typical displacement wins over the spread, and most of 
the particles come out of the box through its one (right) 



E. Exit Time 

The spiting probability is a measure of directionality of 
the first escape from a box of a given width. The quantity 
which is very closely related to the spiting probability is 
exit time, i.e. the time which is needed to escape from 
the same box. Here we study escape time, r, which is the 
first escape time from the box through one of its edges, 
i.e. we do not distinguish between borders of the box. 

The same information which is included in the escape 
time distribution is included in the survival probability. 
Survival probability, S(t) = 1 — F(t), is the probability of 
finding a particle in the box of given width at time t, i.e. 
it is the ratio of particles which at time t are still in the 
box. The survival probability S(t) is exponential function 
of time [la. [20. [361] . In the Fig. \TU\ survival probabilities 
for L = 10 are presented with a — I. Results for other 
values of the box width L and scale parameter a also 
show exponential behavior. 




300 



FIG. 10: (Color online) Survival probability, S{t) = 1-F(t), 
i.e. the probability of finding a particle in the box of the 
width L = 10 for a — 1. Survival probabilities are of the 
exponential type regardless of the value of L and a. 



From first escape time distribution it is possible to cal- 
culate the mean exit time. For a free particle starting its 
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FIG. 11: (Color online) Mean exit time, (r), i.e. the average 
of the first escape time from the box as a function of the 
box half- width for o = 1. Initially a particle is located in 
the middle of the box of the half-width L/2, see Eq. ([2]) and 
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FIG. 12: Value of the exponent b and the prefactor a in Eq. Q 
as a function of the stability index a. Lines present theoretical 
values of the exponent and the prefactor characterizing mean 
escape time for a free particle, i.e. a and 1/[F(q + l)o" Q ] 
respectively. 



motion in the middle of the interval of width L the mean 
exit time is [13 



(r(x(0) = 0)) = 



(L/2Y 



a a T(a + l) 



(8) 



From Eq. |(5J| it implies that (r) oc L a . The same kind 
of behavior can be observed for the Levy ratchet. For 
the large value of the scale parameter a, the presence of 
the potential does not affect the dynamics of the particle. 
Consequently, both the prefactor and exponent are like 
in the free case, see Eq. ((5J. The decrease of the noise 
intensity affects L a scaling and the prefactor in Eq. JSJ), 
especially in situations when the Levy noise is close to 
the Gaussian one and small values of the scale parameter 
a are considered. Nevertheless, for small to moderate 
values of the stability index a, with all considered values 
of noise intensity a, the first escape time is like in the 
free case. Mean exit time as a function of the box half- 
width with u = 1 is depicted in Fig. [TTJ Fig. [12] tests 
applicability of Eq. © . 



IV. SUMMARY AND CONCLUSIONS 

Alltogether, we have numerically documented, that the 
motion of the overdamped particle subjected to a sym- 
metric Levy noise in a potential with a broken spatial 
symmetry leads to occurrence of directionality of mo- 
tion. The preferred direction of motion results as a con- 
sequence of an interplay between non equilibrium charac- 
ter of the underlying stable noise and the broken spatial 
symmetry of the static potential. However, the charac- 
teristics of the particle's displacement, like mean veloc- 
ity and the diffusion coefficient, describing the particle's 
dispersion, get inappropriate when white Levy noises are 



considered. Diverging variance and possibly non-existing 
mean make it necessary to introduce other characteris- 
tics of the directed motion. Consequently, the conducted 
research is based on robust measures characterizing the 
collective motion of Levy-Brownian particles, like group 
velocity, interquantile width of the displacements' distri- 
bution and splitting probabilities. 

In the system under the study, the initially sharp distri- 
bution of particles' positions broadens in course of time 
leading to symmetric probability density of the power- 
law type. The median of the distribution, which charac- 
terizes the group velocity of the particle packet, moves 
linearly with time. The fastest motion of median is ob- 
served for small values of the stability exponent a. The 
coherence of the motion can be characterized by the in- 
terquantile width of the distribution of the particles' po- 
sition, whose growth with time was found to follows the 
1m if) — Qi-mif) °t t 1 / pattern. For a > 1 the width 
of the distribution grows therefore slower than its max- 
imum moves, contrary to the situation when < a < 1 
and the opposite is the case. This fact is mirrored in the 
behavior of splitting probabilities for leaving an interval 
of a given length. 

Eq. (TTJ describes the dynamics of the Levy ratchet for 
any value of the noise parameters. Here, we restricted our 
studies to symmetric noises only. In such a case the only 
relevant noise parameter is the stability index a which 
characterizes the asymptotic power law behavior of ad- 
ditive noise pulses. The stability index controls the level 
of non-equilibrity of the noise. On the one hand, for 
values of a smaller than 2, the Levy noise serves as a 
prototypical perturbation source that drives the system 
out of equilibrium. On the other hand, for a = 2, the 
Levy noise becomes Gaussian noise and models thermal 
fluctuations in the vicinity of the equilibrium state. The 
rectification, ratcheting effect is visible only for values of 
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a < 2. The increase of the stability index a decreases 
the ratcheting effect. Finally for a = 2, similarly to the 
original ratchet and pawl device [J], the absence of the 
average particles' current is a simple consequence of the 
second law of thermodynamics: despite the broken spa- 
tial symmetry assured by the form of the potential, no 
systematic preferential motion of the random dynamics 
can be detected and the ratcheting effect disappears. 

The emergence of the rectification in the "minimal 
Levy ratchet" (Eq. (1)) follows as a result of the non- 
equilibrium type of the stochastic driving and broken 
spatial symmetry of the potential. On the one hand, 
the ratcheting effect is observable for symmetric noises, 
therefore it is induced by the potential. On the other 
hand, in the minimal Levy ratchet scenario the presence 
of the potential is not always explicit: For a free Levy 
particle the probability density of particle positions is 
given by a Levy distribution with growing width and the 
same tails characteristics as the underlying noise. The 
same effect is observed here for the Levy ratchet with 
a < 1.3 when the tail asymptotic of the PDF for a parti- 
cle position becomes characterized by the same exponent 
like the underlying noise. For a ^ 1.3 these exponents 
start to differ suggesting that only for as excessing a cer- 
tain threshold value the particle starts to feel the poten- 
tial. On the contrary, closer examination of the median 
reveals periodic structure of the median position with the 
same space period as the one characterizing the poten- 
tial. Finally, for all values of the stability index a with 
sufficiently large noise intensity a, values of the mean exit 
time scale as in the free particle case. The decrease of 
the scale parameter a makes the system more sensitive 
to the presence of a potential. 

For a Levy ratchet, symmetric stable noise together 



with periodic potential with broken spatial symmetry 
provide a minimal setup for the occurrence of directed 
current. Contrary to Gaussian noises, stable Levy ones 
can be intrinsically asymmetric; the asymmetry of the 
noise introduce preferred direction of the noise pulses 
might affect the current. On the other hand, the asym- 
metry of noise breaks symmetry of the system and intro- 
duces a component which is of 'non-zero average type.' 
The asymmetric Levy noise can be still incorporated to 
the Levy ratchet in the 'zero-average manner.' This can 
be reached by making asymmetry parameter a periodic 
function of time [131 ] . The time dependence of asymme- 
try parameter will slightly modify the type of the Levy 
ratchet making it closer to the ones of either 'thermal' 
or 'tilting' type. Considering such situations might be 
an interesting further step pursuing our line of investiga- 
tions. 
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